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Abstract 

The liquid-gas phase transition for homogeneous symmetric nuclear matter 
is studied in the mean-field approximation. Critical properties are computed 
using a comprehensive group of Skyrme and Gogny forces in an effort to elu- 
cidate the effective interaction dependence of the results. Analytical models 
for the thermodynamical and critical properties are discussed and compared 
to an extensive set of mean-field data. In agreement with these models, a 
tight correlation is found between the fiashing and the critical points. Accu- 
rate predictions for the critical temperature, based on saturation properties, 
can only be obtained after the density dependence of the effective mass is 
properly taken into account. While the thermodynamical properties coming 
from different mean-fields do not follow a law of corresponding states, the 
critical exponents for all the mean-fields have been found to be the same. 
Their values coincide with those predicted by the Landau mean-field theory 
of critical phenomena. 
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1. Introduction 

The liquid-gas phase transition is the canonical example of a phase tran- 
sition in statistical physics At a very basic level, the existence of the 
transition can be directly linked to the structure of the two-body interaction 



Email address: a.rios@surrey.ac.uk (Arnau Rios) 



Preprint submitted to Nuclear Physics A 



July 23, 2010 



[2|. The inter-atomic potential is repulsive at short relative distances and 
attractive at intermediate and long distances. At low temperatures, atoms 
tend to sit at a preferred relative distance from each other, which favors the 
formation of a liquid phase. When temperature is increased, however, the 
thermal motion of the particles is activated and the repulsive, short-range 
structure of the force is probed. The formation of a gas phase is therefore 
facilitated until, above a certain critical temperature, T^, the system com- 
pletely evaporates 

In principle, the classical picture of the liquid-gas phase transition also 
applies to nuclear systems, for which the interaction in the central channel 
has a similar structure 0, 0]. Critical properties might therefore be able 
to provide information on the underlying nucleon-nucleon (NN) interaction. 
In practice, however, the existence of this transition is difficult to establish 
unambiguously ^]. First of all, nuclei are finite systems and the empirical 
signatures of the phase transition might be washed out by their finite size 
effects joj. Moreover, it is difficult, if not impossible, to reconstruct the dy- 
namics of hot nuclei from the available experimental measurements and thus 
the conditions at which these reactions happen cannot be easily determined 
0]. An example of these difficulties is provided by the plateau in the nu- 
clear caloric curve jsf. This phenomenon, observed by different experimental 
groups with various reactions in a relatively wide collision energy range jsj, 
is often taken as a faithful demonstration of the occurrence of the liquid-gas 
phase transition jof. Nevertheless, alternative explanations for the appear- 
ance of such a plateau, that do not resort to any phase transition, can also 



be found in the literature [10|, 111 



Nowadays, experimental efforts are focusing on the study of hot, isospin 



asymmetric nuclear systems |12l . Il3l | . Phenomena such as isoscaling have fu- 
eled the search for connections between the critical point of asymmetric nuclei 
and the isospin dependence of the equation of state. Theoretical analysis of 
these reactions could provide additional constraints to the isospin dependence 
of nuclear forces, besides those given by giant resonances [ij] or transport 



models [15|. An effort in this direction, linking the symmetry properties to 



isoscaling data with a Skyrme-type input for statistical multi-fragmentation 



models, is already underway [16|. Statistical multi-fragmentation models are 
somewhat phenomenologicaL but provide excellent descriptions of several 
observables, such as yields [7(]. 

Less connected to reactions themselves but closer to the spirit of nuclear 
structure studies, there have been several attempts to generalize zero tem- 
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perature calculations to the finite temperature domain. Bulk properties of 
excited nuclear matter can be computed in such approaches, aiming to de- 
scribe the freeze-out conditions and the creation of the primary fragments in 
nuclear collisions. Liquid drop formulae have been used extensively in this 
context 171, ll8|] , in spite of the uncertainties associated with the phenomeno- 
logical finite temperature dependence of their surface and Coulomb terms 



19[. The calculations for isospin asymmetric drops are among the few that 



are able to predict isotopic shifts of the critical properties 20|, |2l|, which are 



of interest given the available experimental data [13 



A more controlled and consistent approach to finite temperature nuclear 
pro perties can be obtained within the self-consistent mean-field framework 
|22| . Mean-field calculations at finite temperature include asymmetry, surface 
and Coulomb effects at once, with no need of additional phenomenological 
parameterizations 23|]. They are based on sound theoretical grounds (sta- 



tistical mechanics together with quantum many-body theory) [2J] and yield 
predictions that are equally valid at zero and finite temperature. Moreover, 
at a formal level, the self-consistent mean- field approach is closely connected 
to density functional theory 22|. The equations and calculations presented 



in the following can also be obtained within the finite temperature extension 



of density functional theory [25|, |26|. The latter theory is among the few 
available frameworks that can provide a controlled methodology to relate 
zero- and finite-temperature predictions, both for nuclei and nuclear matter. 
So far, calculations for hot nuclei have been generally performed on selected 
spherical isotopes, mostly aimed at understanding the structural changes in- 



duced by temperature |27|, |28|, |29| . In principle, finite temperature mean-field 
calculations could also be used to find the limiting temperatures for selected 
isotopic chains. The necessary extension to deformed isotopes should be fea- 
sible with present computational capabilities. Yet, in spite of their great 
interest, systematic studies of this type for a relevant range of isotopes are 
still lacking. 

Before such a program is implemented, however, the model dependence 
of the calculations should be assessed. Key among the possible sources of 
uncertainty is the dependence associated with the underlying NN effective 
interaction. It is well-known that different effective interactions predict dif- 
ferent bulk properties for nuclear matter 30|. This statement also applies 
at finite temperature and, as a result, different liquid-gas critical points are 
found for different forces. In theory, a reliable determination of the critical 
point could be helpful in constraining the properties of the effective interac- 
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tion [3l|. In this contribution, I will discuss the case of symmetric nuclear 
matter, where this connection should be more straightforward due to the 
absence of symmetry. Coulomb and surface terms. 

Because of its simplified nature, the mean-field picture is presumably 
unable to describe correctly the whole phase diagram of nuclei or nuclear 
matter. The level density lacks the effect of complex many-body excitations 
and therefore underestimates that of nuclear systems 32]. Even in infinite 



matter within a fully quantal description, the effect of beyond mean-field cor- 
relations and clusterization on the density of states and the thermodynamics 
of the system is difficult to quantify. As for finite nuclei, clusterization is 
completely absent in the mean-field picture, but it is known to play a major 



role in hot nuclear systems [33[. In general, the mean- field approximation 
provides a poor description of the critical properties of finite, charged sys- 
tems Related to that, information gathered in time- dependent studies 
indicates that the mean-field description gives a particularly poor description 



of the spinodal zone |3J, |35[. The fragmentation mechanism of nuclei is sen- 
sitive to many-body correlations and dynamical time-scales in the reaction 
dynamics, and the mean-field picture fails in describing such effects. 

Consequently, the mean-field description of nuclear matter should not 
be taken as meaningful description of the liquid-gas transition in nuclear 
systems. As a matter of fact, the study presented here does not try to 
provide a quantitatively accurate picture for such critical properties. Instead, 
the aim of this study is to investigate the effective interaction dependence 
of the liquid-gas transition properties. As such, this can be taken as a first 
methodological contribution to test what information, if any, can be obtained 
from such analysis. Homogeneous nuclear matter within a mean-field picture 
provides the simplest testing ground for such an analysis due also to the 
large number of available effective interactions. Other theoretical approaches, 
which, for instance, provide a more realistic description of the spinodal zone, 
will have an underlying dependence on the effective interaction and similar 
analysis can be carried out for these cases. 

In addition to providing a methodological test, the study of the critical 
properties of homogeneous matter is also of theoretical interest in connec- 
tion with recent developments. First of all, the liquid-gas phase transition 
can nowadays be computed from many-body approaches that go beyond the 
mean-field. Starting from realistic phase-shift equivalent NN interactions and 
using finite temperature many-body techniques to find the thermodynami- 



cal properties, the liquid-gas phase transition can be constructed [36|. Such 
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ab initio calculations provide a relatively wide range of predictions for the 
critical properties dependin g on the underlying interaction and many-body 



approximation that is used [37|, |38| . It is therefore important to investigate 
whether a similar wide range of values is obtained within mean-field theory 
and whether such an interaction-dependence of the results can be eliminated. 

Secondly, the phase transition for homogeneous matter has important 
connections to astrophysics, particularly to the onset of instabilities in the 



isospin asymmetric nuclear medium |39l. |40| . In this context, one would like to 



know whether the onset of the instabilities is strongly dependent on the prop- 
erties of the underlying interaction. Astrophysical observations could then, 
in principle, provide constraints for nuclear effective interactions. Moreover, 
with over a hundred Skyrme force parameterizations at hand, it is nowadays 
possible to carry out "statistical" calculations with a comprehensive set of 



different mean- fields [41[. If one supposes that all these forces are equiva- 
lent, their predictions should span a physically relevant parameter space of 
observables within the mean-field approximation. Focusing on the liquid-gas 
phase transition, the predictions obtained for different interactions should 
provide a set of vahd theoretical critical points. Similar studies, relating 
the liquid-vapor coexistence curve to the underlying interaction, have been 



performed in the context of chemical physics j42 



By computing the liquid-gas critical properties with large sets of mean- 
fields, one might also be able to identify behaviors which are independent of 
effective interactions. Examples of such correlations will be given in Section 
m In addition, a comprehensive study of the critical points predicted by 
forces with different saturation properties can assess in a systematic way the 
long-standing question of whether there is a connection between the criti- 
cal and the saturation points. A link between these two physically relevant 
points might yield relevant information about the temperature dependence 
of the equation of state. Moreover, extending these results to the case of 
isospin asymmetric matter could provide a connection with recent studies 
on correlations between zero-temperature, isospin-asymmetric nuclear prop- 



erties 



41| . Finally, in the context of critical exponents, one would expect 
these to be independent of the underlying mean-field. This can be tested nu- 
merically with self-consistent calculations of the liquid-gas phase transition, 
providing a nice example of how physics at the critical point is dominated 
by fluctuations. 

This paper is structured as follows. Section |2] provides a brief summary of 
the formalism associated with the finite temperature self-consistent Hartree- 
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Fock. The properties of the hquid-gas phase transition and the spinodal and 
coexistence curves are also introduced. In Section [31 analytical models that 
provide a useful qualitative understanding for the liquid-gas phase transition 
are discussed. The relation between the saturation, flashing and critical 
points of nuclear matter is explored in Section HI Section [5] is devoted to 
analyze the law of corresponding states and the critical exponents of the 
phase transition for self-consistent mean-field calculations. The conclusions 
and outlook for future work are outlined in Section |6l 



2. Self-consistent mean- field calculations of hot nuclear matter 



The self-consistent Hartree-Fock method for nuclear matter calculations 
offers a theoretically consistent model for the description of hot nuclear mat- 
Physically, this formalism stems out from a minimization of the 
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ter 

free energy of the system subject to the condition that the density matrix is 



of the one-body type j23|. Two- and higher-body correlations are therefore 
neglected on the basis that their contributions are captured by the parame- 
terization of the effective interaction. Alternatively, one can cast the problem 
in terms of density functional theory, for which an energy density functional 
is guessed (based on some knowledge of the system) and a corresponding local 
density is found via a minimization of the energy (or, at finite temperature, 
the grand-potential |25l. |26|). 

In the following, results are presented for the nuclear matter liquid-gas 
phase transition with two different types of widely used effective interactions. 
The Skyrme force is a density-dependent, zero-range parameterization: 



XiPs 



•(—2 . . — >2- 

k 6{'^) + 6{l^)k 



+ t2{l + X2Ps) ^ ■S{l^)k 



(1) 



where is the relative distance between two interacting nucleons and k 

2)/22 is their relative momentum acting on the right 
conjugate acting on the left, while Ps = {1 + ■ a^2)/2 represents the spin 



k is its 



exchange operator and p corresponds to the total nucleonic density |43| . The 



spin-orbit term is disregarded, since it does not contribute in homogeneous 
systems. Over a hundred parameterizations of this force exist in the literature 
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30|. Results for 39 forces will be presented. 27 of these passed a series of 



criteria imposed in Ref. [30|: they offer a good description of both nuclear 
and neutron matter. In addition, to provide checks with previously published 
studies, the SIII, SkM*, RATP and SkS3 parameterizations j4J], together 



with 2 forces of Ref. |l8| (forces SI and SJl) and Ref.[5J] (forces ZRl and 
ZR2), have also been considered. Finally, more recent mean-field models, like 
the LNS parameterization jisl, fitted to microscopic Brueckner-Hartree-Fock 



calculations, as well as the latest force of Goriely et al. [46|, which provides 
an excellent description of finite nuclei properties, have been included in the 
discussion. 

The Gogny interaction has also been used extensively for mean-field and 
pairing calculations of nuclei and nuclear matter. In addition to a density 
dependent zero-range term, it includes a set of finite-range terms: 



1=1 

+ to(l + a;oP.)p"5(^). (2) 



as well as a more complex spin-isospin structure (as seen by the additional 
isospin exchange operator, Pt). There are far fewer Gogny forces than Skyrmes 
in the literature. The following parameterizations will be considered: Dl (47| . 
DIS [48], DIP (with its additional zero-range term) [i^ and the very recent 



DIM force 50|. The D250-D300 parameterizations of Ref. 5l|] will also be 
included in this analysis, since they might be useful in highlighting a relation 
between the compressibility and the critical properties of Gogny forces. 

In contrast to the zero temperature case, for which the equations for 
the microscopic and macroscopic properties of nuclear matter can be cast 
analytically for both effective interactions, the self-consistent Hartree-Fock 
approximation needs to be implemented numerically at finite temperature. 
The formalism presented in the following is grand-canonical and the external 
fixed variables should be the temperature, T, and the chemical potential, fi. 
In the thermodynamic limit, however, ensemble equivalence is guaranteed 
and one can work equivalently in a canonical picture, where the average 
density of the system, p, instead of /x, is considered as an external variable. 
The energy density in the Hartree-Fock approximation is given by the sum of 
the kinetic term plus a potential term, obtained from the expectation value 
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of the effective interaction between two plane waves of momentum ki and ^2 '■ 



(p,T) = u / 77T-T^/(A;)- 



+ y / ^|^/(fci)/(^2)(fc";fc;i^;v^ifc";fc2)a, (3) 

where z/ = 4 is the spin-isospin degeneracy of the system and the properly 
anti-symmetrized matrix elements have been averaged over spin and isospin. 
The momentum distribution, f{k), describes the thermal population of mo- 
mentum states according to a Fermi-Dirac distribution, 



fik) 



1 + e 



'!^+u{k)-^;)/T 



(4) 



and depends on both the chemical potential, fi, and the mean-field, U{k). 
The chemical potential can be found by inverting the relation: 

which yields the normalization of the distribution to the total density of the 
system. The mean-field is given by: 

U{k) = uj ^J{k'){kk'\VMN\kk')a, (6) 

and also depends on the momentum distribution, f{k). As a result, Eqs. ([5]) 
and ([6]) are coupled and need to be solved simultaneously, which gives rise to 
the self-consistent Hartree-Fock solution. Once the momentum distribution 
is set, it can be used to compute bulk thermodynamical properties, such as 
the entropy density: 

§(p,T) = z. j ^[/(fc)ln/(A:) + (l-/(fc))ln(l-/(A:))]. (7) 

From this and the energy of Eq. (jS]), the free energy is obtained, F = E — TS, 
and access to the pressure is gained by taking a derivative of the free energy 
per particle, p = P^^f^- To avoid the latter differentiation, one usually 
replaces the derivative with the chemical potential, so that p = p{fi—^). This 
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replacement relies on the thermodynamical consistency of the theory, which 
is only achieved if the rearrangement terms of the mean-field, originated by 
its density dependence, are properly considered in the inversion of Eq. ([5]) 
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The solution of the Hartree-Fock equations is particularly easy for the 
Skyrme mean-field, due to its simple structure. Most importantly, the mean- 
field is independent of temperature and only depends quadratically on mo- 
mentum. This allows for a separation of the type U{k) = U + ^^^Ap. The 
first term, including the rearrangement contribution, amounts to: 

U=ltop + ^{a + 2)tsp^+\ (8) 

while the momentum dependent term can be easily incorporated in the defi- 
nition of the effective mass: 

m 1+ 2m^(3ti + 5t2 + 4^2X2) p 1 + Ap' 

Note that, for Skyrme forces, the mean-field is independent of f{k), U can be 
absorbed into the chemical potential and the effective mass contribution can 
be included in the kinetic term. As a consequence, the numerical solution of 
Eq. ([5]) is substantially simplified. In particular, even though at the diagram- 
matic level the calculation is self-consistent, in the numerical implementation 
there is no need for a self-consistency loop. 

The determination of the mean-field is more difficult with the Gogny 
force. In this case, due to the finite range of the interaction, the single- 
particle potential is momentum, density and temperature dependent: 

3/2 2 

U{k, p, T) = ^ ^^ [4^^ + - '^H, - M,] 

i=l 

2 

- n'/' Yl + - 2^. - g,{k, p, T) 
1=1 

+ ^to{a + 2)p"+\ (10) 
The exchange integral, 

9.{k,p,T) = . J ^,e-tr'-'y{k'), (11) 
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Figure 1: Pressure isotherms for different temperatures: T = (solid line), T — Tf 
(dashed Une), T = (dash-dotted Hne) and T = 18 MeV (dash-double-dotted line). The 
dotted lines represent the corresponding spinodal and coexistence regions. The left panel 
has been obtained with the BSkl7 interaction, while the right panel shows results for the 
Gogny DIM force. 



is responsible for the non-trivial temperature and momentum dependences. 
The Gogny mean-field cannot be separated cleanly into an effective mass 
and a momentum independent term. This separation might eventually be 
performed as an additional approximation and, on the whole, should pro- 
vide reasonable results |52i]. However, full mean- field calculations are prefer- 
able for accurate calculations of liquid-gas coexistence curves. Because of 
the implicit dependence of the exchange integral on the mean-field, a self- 
consistency loop at the computational level, absent in Skyrme calculations, 
is needed. 

The density and temperature dependence of the thermodynamical prop- 
erties obtained within the Hartree-Fock approximation are non-trivial. In 
particular, they cannot be parametrized in analytical form except for extreme 
regimes of degeneracy. The mean-field solution describes simultaneously two 
different physical phases: a gas at low densities and a liquid at densities close 



to saturation [28|. This behavior is more easily visualized in the pressure 
isotherms of nuclear matter, as show in Fig. [H Left and right panels corre- 
spond to results obtained with typical Skyrme (BSkl7) and Gogny (DIM) 
forces, respectively. Four characteristic temperatures have been chosen for 
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illustrative purposes. 

At T = (solid lines) , the pressure is a decreasing function of the density 
at low densities until a minimum at ~ Po/2 (where po is the saturation 
density) is reached. This regime, for which ^ < 0, corresponds to the 
spinodal region, where the system is mechanically unstable. As temperature 
increases, the spinodal region, which is shown in light dotted lines in Fig. [H 
shrinks. At low densities, a mechanically stable gas phase forms and coexists 
with the higher density liquid phase. The spinodal region is associated with 
negative pressures at low temperatures. This picture changes as temperature 
increases and the thermal contribution tends to push the pressure to largest 
values. The isotherm at which the pressure becomes a completely positive 
function defines the so-called flashing temperature. More specifically, the 
flashing point satisfies the two simultaneous constraints: 

(Note that, for nuclear matter, the point at T = and p = also satisfies 
this condition). In Fig. [H the flashing point corresponds to the crossing of 
the spinodal with the p = line. For an isolated system, without an exter- 
nal gas to stabilize it, this would correspond to the maximum temperature 
at which the system could still be self-bound The definition of the flashing 
temperature, T/, is therefore somewhat similar to that of the limiting tem- 
perature of a nucleus, i.e. the maximum temperature a nucleus can withstand 
before evaporating into a gas of nucleons. As a matter of fact, Tf might be 
considered as the extrapolation of the limiting temperature to the ^ ^ oo 
case. 

For the infinite system, however, the system can exist above the flashing 
temperature because the gas and the liquid phase exert the same pressure 
on each other. According to the Maxwell criterion, the condition that the 
gas and liquid phases have the same pressure, together with the requirement 
that their chemical potentials are the same, set the coexistence points [H, lij. 
The coexistence region is shown with a dark dotted line in Fig. [H Note that, 
at T = 0, there is a coexistence between a zero-density gas and a liquid at 
saturation density. Increasing the temperature results into larger gas and 
smaller liquid densities. Eventually, the two densities meet at the critical 
point, which occurs at the critical temperature, Tg. The critical isotherm is 
shown with a dashed-dotted line in both panels of Fig. [H The critical point 
is found at the endpoint of both the spinodal and the coexistence lines, and 
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therefore can be found by requiring: 



Above the critical temperature, the system will only exist in the gas phase 
and the pressure isotherms become a monotonically increasing function of 
density. An example of this behavior is given by the dashed-double-dotted 
lines, which have been computed at T = 18 MeV for the two panels. 

For optimal liquid-gas calculations, it is important to compute Fermi in- 
tegrals [integrals over the Fermi-Dirac distributions such as those in Eqs. ([3]), 
(I5l), (E]) and with great accuracy and, simultaneously, short computa- 



tional times. Numerical techniques based on Ref. [53| have been implemented 
to this end. For consistency, the critical point has been determined by three 
different numerical methods. The first two estimates have been obtained by 
finding the endpoint of the coexistence and the spinodal curves, respectively. 
These curves are computed separately using Newton-Raphson routines and 
extrapolations are performed to determine their maximum accurately. In 
addition, the last temperature and density where the first and second deriva- 
tives of the pressure show a node has been found using a bisection method. 
The three different methods coincide in the determination of the critical tem- 
perature to the second decimal point for all the mean-fields considered here. 
For the forces where it was possible to do so, the results presented have 
been checked against those previously published in the hterature. To my 
knowledge, this is the first time that a consistent and comprehensive set of 
mean-fields is used to find correlations between physical parameters in the 
liquid-gas phase transition. 



3. Models for the phase transition 

Finite temperature nuclear studies provide new means of understanding 
the properties of the dense nuclear medium. With the suitable theoretical 
framework, the information gathered from hot nuclei could be used to con- 
straint the nuclear mean-field. In particular, two points of the phase diagram 
of nuclear matter stand out for their relevance. On the one hand, the sat- 
uration point defines the properties of cold nuclear matter and is directly 
related to nuclear structure data. On the other hand, the liquid-gas critical 
point can only be accessed if excited matter is produced in a reaction. Can 
one find a direct connection between the critical and saturation points? In 
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Figure 2: Left panel: saturation curve of nuclear matter for four different interactions. 
Right panel: liquid-gas coexistence curves for the same interactions. 



principle, the partition function of an interacting system is not constrained 
by its zero temperature properties. For nuclear matter within the mean-field 
picture, however, the effective interaction determines the energy and also 
the entropy [via the single-particle potentials in Eq. ([7])] and perhaps this 
indirect connection might be strong enough to provide a correlation between 
such properties. 

At first sight, however, such a connection is not straightforward. An 
example is provided by Fig. [21 where the saturation curve (left panel) and 
the liquid-gas coexistence curve (right panel) are shown for four different 
mean-fields. The Skyrme forces ZRl and SLy9 have the highest and lowest 
critical temperatures, respectively, of the interactions considered here. DIM 
and BSkl7 should be taken as representatives of the "average" behavior of 
Gogny and Skyrme forces. Even though these four interactions have rela- 
tively similar saturation points, with saturation energies deviating by less 
than 1 MeV from one another, the corresponding critical temperatures span 
a wide range of values, from 14 to 23 MeV. This suggests that no connection 
exists between the saturation and the critical point. 

In contrast to this naive analysis, data concerning nuclear caloric curves 
has been used to provide constraints on the properties of nuclear matter at the 
saturation point [31]. The knowledge gathered in this direction, however, has 
been generally based on models that connect saturation and critical prop- 
erties by assuming analytical density and temperature dependences of the 
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equation of state. By construction, the thermodynamical properties within 
these models have underlying temperature and density dependences that 
might or might not be good approximations of the full nuclear system. Gen- 
erally, such models are based on expansions in particular degeneracy regimes. 
The degeneracy parameter, d = j^, where T is the temperature and €p, the 
Fermi energy, measures the relative importance of thermal effects against 
quantum fluctuations. The thermodynamical properties in the Hartree-Fock 
approximation can be found analytically at the extremes of degeneracy. On 
the one hand, in the degenerate regime {d << 1), i.e. at relatively low tem- 
peratures or high densities, the Sommerfeld expansion provides a systematic 
expansion of the thermodynamical properties in terms of c? j2|. This yields 
the well-known quadratic (linear) temperature dependence for the energy 
(entropy). On the other hand, in the semiclassical limit [d » 1), i.e. for 
very dilute and hot systems, the fugacity expansion provides a connection 
with classical statistical mechanics 

Notions of both the degenerate and the semiclassical limits have been 
used in the literature to model the behavior of the hot equation of state 



54 |55[. Unfortunately, the actual nuclear liquid-gas phase transition takes 
place at T ~ 15 MeV and p ~ 0.08 fm~^, so that c? ~ ^ and none of 
the above expansions can be applied with extreme confidence. Nowadays, 
full self-consistent calculations can be implemented numerically and very 
accurate results are found in extremely short computational times (seconds 
for the case of Skyrme forces) . Consequently, there is no need to rely on such 
expansions other than to get a physical insight on the properties of liquid-gas 
transition. In the following, I shall discuss some of these models and analyze 
their validity for the liquid-gas phase transition. In particular, I will look at 
the relation between saturation and critical parameters. Incidentally, these 
models will also provide insight into a direct connection between the flashing 
and the critical points, as will be discussed in Section HI 

3.1. Jaqaman model 

A very simple model that already encodes some of the physically relevant 
information associated with the nuclear liquid-gas transition is provided by 
the lowest order fugacity expansion with m* = m. This model (and higher 
order contributions in the expansion) were discussed by Jaqaman et al. in 



Ref. 5J]. Within this approach, the temperature dependence of the pres- 



sure is taken to be that of a classical gas and the density dependence is 
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supplemented with the zero-range terms of a Skyrme force: 



p = pT + iop' + ^t3(« + l)p"+' 
8 lo 



(14) 



This expression neglects any quantum fluctuation (in the sense that the tem- 
perature dependence is purely classical) and, strictly speaking, should only 
be valid in the semiclassical regime. The flashing point can be obtained 
analytically from Eq. (fT2|) : 



Pf 
Tf 



_6 to 

L(a + 1)2^3 

1 a 



l/a 



(15) 
(16) 



while the critical point, deduced from Eq. ( |T3|) . reduces to: 



Pc 
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_(a + l)2(a + 2)t3 
3 a 

toPc , 



l/a 



_ a + 1 



(17) 
(18) 

(19) 



This expression gives an important insight into the relation between the 
critical point and the underlying equation of state. The critical density is 
determined by the ratio of to and ts, so it will be sensitive to the interplay 
of the attractive and the repulsive components of the force. Moreover, pc 
depends non-analytically on a. The latter point stresses the importance of 
the density dependence of the force in determining the critical point. In 
particular, the compressibihty, Kq, is to a great extent fixed by a and, thus. 



it will play a role in the determination of the critical point [3l|, |5l|]- The so 



called critical compressibility factor is given by the dimensionless parameter 
7c = is useful in characterizing the critical point [ij. For the 

van der Waals equation of state, 7c = | = 0.375 while, for the Jaqaman 
model, 7c = 2{a+2) • -^^^ Skyrme forces, which usually have a = — 1, this 
corresponds to lower critical compressibility factors than the van der Waals 
prediction, between | < 7c < |- 
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Within the Jaqaman model, the critical and the flashing points are closely 
connected. Their ratios are purely determined by the exponent of the density 
dependence: 



pc 

Pf 



a + 2 



l/a 

(20) 



^ = 2^. (21) 
^/ Pf 

For Skyrme mean-flelds, this results in a very limited range of allowed ra- 
tios: ^ < ^ < |. Accordingly, the quotient between critical and flashing 

temperatures is twice as large, ^ ~ 1.3. 

The Jaqaman model provides a very poor description of the liquid branch. 
The latter involves densities close to saturation and therefore sits closer to 
the degenerate regime. Speciflcally, the linear temperature dependence of 
the pressure is only valid at extremely low densities and therefore it is too 
simple to explain accurately the whole region involved in the transition. This 
simplistic parameterization of the finite temperature properties leads to large 
differences between the critical temperature of the full self-consistent calcu- 
lation and those of the respective Jaqaman models (see Fig. |5]). 

3.2. Kapusta model 

The expansion near the completely degenerate case can also be used to 
establish the temperature dependence of the equation of state. The Kapusta 
model, originally derived in Ref. ^5|, relies on the Sommerfeld expansion 
together with the knowledge of the compressibility to build an approximation 
for the pressure of the system: 

P=^4(P-Po) + ^&y/='T^ (22) 

Ingredients of this model include the numerical constant b = ( ) and 

the effective mass at saturation density, rriQ = m*{po). The density depen- 
dence of the latter is neglected, since ml is basically introduced to tune the 
thermal density of states close to saturation. Again, analytical formulae can 
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be found for both the flashing: 




1/3 



and the critical points: 



pc = —po , 
55/6 



Pc 



213/631/35 

25 




1/3 



324 



Kopc 



(23) 
(24) 

(25) 
(26) 
(27) 



Plugging in typical values of saturation parameters (po = 0.16 fm ^, Kq ~ 
230 MeV and mg = 0.7m), the critical point is predicted to be at pc = 
0.066 fm-3, Tc = 20.7 MeV and pc = 1.18 MeVfm-^. While pc and are 
rather typical and comparable to mean-field values (see Table [1] below) , Pc 
is somewhat large. This suggests that the pressure isotherms in this model 
overestimate the mean- field ones. As a result, the critical compressibility 
factor becomes extremely large, 7c = 0.89, in contrast to the van der Waals 
and mean-field values. 

Within the Kapusta model, both pf and pc are directly proportional to 
the saturation density. This suggests that a clean connection between critical 
and saturation densities might be established. Mean-field calculations do not 
support these findings, as will be discussed later on (see Fig. [2]). Just as in the 
previous model, however, the critical and flashing points are closely related 
to each other: 



pc 
Pf 



2 

3 



T 



/ 



16 y/^ 
3 / 



(28) 
(29) 



The ratio of densities coincides with the upper value predicted by the Jaqa- 
man model. The flashing and critical temperatures are functions of the com- 
bination 1] = 



1 /3 

pQ , SO that their ratio is a simple numerical factor. As 
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a matter of fact, ^ = 1.75, which is substantially larger than the value 
predicted by the Jaqaman model. 

The dependence of Tf and Tc on 77 is helpful in understanding the links 
between the saturation and critical points. To begin with, stiffer equations 
of state lead to larger critical temperatures. At T = 0, the pressure in the 
Kapusta model has a minimum value of Pmin = — ^-^oPo at sub-saturation 
densities. High values of the compressibility or larger saturation densities are 
therefore associated with more negative pressures below saturation. Conse- 
quently, higher temperatures are needed before the pressure can become zero 
(flashing point) or its minimum disappears (critical point). This explains why 
Tc and Tf are monotonically increasing functions of Kq and po- Similarly, 
"stronger" temperature dependences (in the sense that the density of states 
is larger due to an increase in mg) can overcome negative pressures easily. 
Larger effective masses are therefore associated with lower critical and flash- 
ing temperatures. 

The Kapusta model is based on a degeneracy expansion which is only 
valid at low temperatures and large densities. As a result, the pressure 
isotherms of Eq. fl22|) are very close to the mean- field ones at T = 0. As 
temperature is increased and the critical point is approached, however, the 
agreement worsens. The origin of these discrepancies can be traced back 
to the simplified temperature dependence, but also to the lack of density 
dependence in the density of states. The former effect can be improved by 
considering higher order terms in the degeneracy expansion, while the latter 
is likely to play an important role for the liquid-gas phase transition, since the 
flashing and critical points are largely determined by the density derivatives 
of p [see Eqs. ( |T2l) and ( |T3|) ]. An extension of the Kapusta model which 
includes a more realistic effective mass dependence is discussed in Section 



3.3. Lattimer-Swesty and Natowitz predictions 

Two popular parameterizations of the critical temperature have been in- 
spired by the Kapusta model relation between Tr and the saturation proper- 
ties. Lattimer and Swesty introduced in Ref. [56| a formula for the critical 
temperature in terms of saturation properties: 



While based on the Kapusta model, this expression has the opposite de- 
pendence on saturation density i.e. larger saturation densities lead to lower 



in 




(30) 



18 



critical temperatures. The value for the numerical parameter, C^s = 0.608 
MeV^/^fm~^, is obtained by fitting the critical temperature of a series of 



simple mean-field models [56 



Later on, Natowitz et al. generalized the previous expression to include 
the effect of the effective mass near saturation, ml: 



The empirical value of the parameter C^v = 0.484 ± 0.074 MeV^/^fm"^ is de- 



termined from an average of several theoretical values [31|. In contrast to the 
Jaqaman or the Kapusta models, which try to describe the thermodynamical 
potential, these two predictions are phenomenological estimates that provide 
a fit to a series of critical point data. The extent to which these predictions 
is satisfactory will therefore depend on the original data used in the fits. 
In the following Section, the accuracy with which these formulae reproduce 
the nuclear matter mean-field results with Skyrme and Gogny forces will be 
discussed. 

3.4- Extended Kapusta model 

The last model to be presented here is an extension of the Kapusta model 
that, to my knowledge, has never been considered. It provides a particularly 
good prediction of the critical temperature in terms of saturation properties 
for mean-field calculations. A common feature of the models presented so far 
is the fact that they completely neglect the density dependence of the effective 
mass. The latter changes from values of around — = 0.7 at saturation to 1 
at p = 0. Since the effective mass sets the density of states at each density, 
this (up to) 30 % effect might be particularly important for the density 
and temperature dependence of the equation of state and, consequently, for 
accurate predictions of critical parameters. Moreover, the effective mass 
(particularly its energy dependence) has been acknowledged as one of the 
key parameters in determining nuclear caloric curves 

To include this effect schematically, one can use as a guiding point the 
density dependence of m*{p) coming from the Skyrme interaction, Eq. (|9]). 
This dependence enters the Kapusta model in the thermal term and modifies 
the expression of the pressure: 

Kop\ m 2 1 + |AP 1/3^2 fno\ 
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Note that this differs exphcitly from Eq. (122 p even at p = po- In the original 
derivation of the Kapusta model, the only density dependence in the term 
comes from the Fermi energy in the denominator. The extended model in- 
cludes more complicated dependences that arise from the density derivatives 
associated with the calculation of the pressure. 

Given this expression, one should go ahead and keep the density depen- 
dence in the solution of the flashing and the critical point equations. Unfor- 
tunately, the latter task is complicated by the fact that analytical solutions 
cannot be found. However, since both the flashing and critical points occur 
well below saturation density, the effective mass is close to the free mass and 
the dimensionless parameter Ap is of order 10 % or so. One can then expand 
all the expressions in terms of this parameter and solve the corresponding 
equations to an accuracy of (9(A^p^). Using this method, the flashing point 
is found to be: 



9 

Pf = Pf 1 - ^Ap/ 
Tf =Tf[l- Apf] , 



(33) 
(34) 



where the tilded parameters correspond to those of the original Kapusta 
m. Similarly, the critical point becomes: 



model with ttIq 



Pc 

T 

Pc 



pc 



9 

1 - — Ap, 
20 ^ 



Tc [1 - Ap,] , 

Pc [1 - Apc] . 



(35) 

(36) 
(37) 



For most Skyrme parameterizations, the effective mass decreases with 
density, so that A > 0. As a consequence, both the critical and flashing den- 
sities are reduced by a very small numerical factor. The effective mass den- 
sity dependence is more important on the flashing and critical temperatures. 
Note that, interestingly enough, and up to first order, both temperatures can 
be written as: 



Tr 



m*iPf) 

m 
m*{pc] 

m 



T 



(38) 
(39) 
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These expressions suggest that the effective mass below saturation is impor- 
tant for a proper determination of the flashing and critical temperatures. 
Generally, for Skyrme forces both effective masses are around ~ 0.85 — 0.9, 
which leads to a sizable effect, particularly if one is looking at accurately 
predicting critical values from saturation properties. The Skyrme force re- 
sults presented in the following have been obtained with Eq. fl5^ . Moreover, 
these expressions have the virtue of being easily generalizable to effective in- 
teractions that have other (presumably more complex) density dependences. 
Assuming that the temperature dependence of the effective mass can be 
neglected, as is generally the case for Gogny forces, one can compute the ef- 
fective mass at the Fermi surface at T = and plug it into these expressions. 
The examples with the Gogny force discussed in the following have been ob- 
tained in this way. Finally, these formulae can also be used in calculations 
performed with sophisticated many-body theories. 

4. Relation between saturation and critical properties 

In this Section, the relation between the saturation and critical parame- 
ters predicted by self-consistent mean-field models will be analyzed by means 
of systematical calculations with an extensive set of effective interactions. 
Critical and flashing points will be deduced and compared to the predictions 
given by the models of the previous section. Thanks to the large number 
of mean-field parameterizations at hand, correlations between different pa- 
rameters will be identified. With this, I will highlight the advantages and 
disadvantages of the different models and will try to draw a coherent picture 
of the liquid-gas transition properties as derived from accurate self-consistent 
calculations. In addition, this analysis is useful in assessing the correctness 
of naive dimensional arguments that have been put forward in the context of 
the liquid-gas phase transition. Note, however, that these correlations might 
not be physically realistic, since the description of the liquid-gas transition 
close to the critical point within a mean-field picture can be put into question 

53- 

The upper panels of Fig. |3] concentrate on the relation between the flash- 
ing and the saturation points. The flashing points have been determined 
for several Skyrme (circles, triangles and squares) and Gogny (diamonds) 
forces. For reasons that will become clear later on, Skyrme interactions with 
no effective mass (^ = 1 or A = 0) are shown with squares. The circles 
correspond to modern Skyrme forces, i.e. those with m* ^ m that passed 
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Figure 3: Left panels: flashing (upper panel) and critical (lower panel) densities versus 
saturation density for different effective interactions. Right panels: flashing (upper panel) 
and critical (lower panel) temperatures versus saturation energy. Skyrme interactions 
with 2i- = 1 are represented with squares. The circles correspond to results obtained with 
modern Skyrme forces that have ^ ^ 1, while older forces are shown with triangles (see 
text for details). The diamonds have been obtained with the Gogny interaction. 



the restrictions of Ref. 30|] or those that have been introduced later on (LNS 
and BSkl7). Triangles correspond to somewhat older forces (SI, SIII, SkM* 
and RATP). The upper left panel shows a scatter plot of the flashing densi- 
ties obtained with finite temperature self-consistent calculations for different 
forces against their respective saturation densities. Any correlation between 
these parameters should appear as a clear trend in these figures. 

In most models, the saturation density is set (or fitted, depending on the 
particular model) to a value between 0.145 and 0.175 fm^^. The range of 
allowed flashing densities is relatively narrower and goes from around 0.08 to 
0.10 fm~^. Actually, all p/'s fall between two "empirical" correlations lines. 
Limiting by below is the relation pf = 0.5pQ, whereas the upper limit is set 
by the Kapusta model prediction, pf = 0.625po- At first sight, however, it 
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is hard to find a clear proportionality, as suggested by such relations. As a 
matter of fact, forces with the same saturation densities lead to fiashing den- 
sities that differ by more than 0.02 fm~^, suggesting that additional physical 
information is needed if both densities are to be connected. 

A somewhat similar situation is found by comparing the flashing temper- 
ature and the saturation energy, cq (upper right panel). Naively, one could 
have argued that the temperature needed to vaporize completely a liquid 
might be proportional to the binding energy of the latter. In contrast to 
this, none of the models presented in the previous Section predict an explicit 
dependence on cq. On top of that, no clear proportionality trend between 
Tf and Cq is observed in Fig. El Actually, for all forces, the flashing temper- 
atures are well below the corresponding saturation energies. Furthermore, 
for mean-flelds that have very similar saturation energies, the corresponding 
flashing energies differ by more than 4 MeV. 

The lower panels of Fig. |3] focus on the connection between the critical 
and the saturation densities. Similarly to what is observed for the flashing 
density, most of the pcS are found in a narrow range of values. A band, 
bound by below by the "phenomenological" law, = |po, and by above by 
the Kapusta model prediction, pc = ^po = OAlGGpo, describes well most 
of the data. Yet, it is hard to extract a clear proportionality law from this 
picture. The lower correlation seems to describe particularly well some of the 
critical-saturation points, but one should not take the correlation pc = |po as 
a faithful prediction for critical densities. Similarly, looking at the lower right 
panel of Fig. |3l some critical temperatures happen to lie close to the Tc = cq 
line, but this law is hardly generalizable to all effective interactions. An even 
more striking fact is the noticeable spread in T^. while all saturation energies 
lie within 2 MeV from each other, critical temperatures span a range of almost 
9 MeV. This indicates that there is no connection between the saturation 
energy and the position of the critical point. It is interesting to note that the 
magnitude of the spread in is somewhat similar to the differences in critical 
temperatures found by using different realistic interactions and many-body 
techniques in Ref. [37 . 

The scatter plots in the upper and lower panels of Fig. Elshow very similar 
structures, with a clustering of Skyrme values around the central region, a 
more dilute set of Gogny forces and some m* = m forces in the right end. It 
appears that subsets of forces with analogous saturation properties behave 
similarly in their prediction for critical and flashing points. Accordingly, 
critical and flashing properties for different forces could be correlated. Such 
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Figure 4: Left panel: critical density versus flashing density for different effective interac- 
tions. Right panel: critical temperature versus flashing temperature for different effective 
interactions. 



correlations are explored in the two panels of Fig. HI In the left panel, critical 
densities obtained with different mean-fields have been plotted against their 
corresponding flashing densities. As suggested by both the Jaqaman and 
the Kapusta models [see Eqs. (120|) and (!28l) ]. a clear linear proportionality is 
found between the two. Fitting this linear relation with a single parameter, 
the law pc = 0.65pf is found, in good agreement with the proportionality 
constants predicted by the Kapusta and Jaqaman models. A more accurate 
fit of the larger densities can be achieved by considering an offset. The 
corresponding linear regression yields the relation = —0.02 + 0.89p/. 

Similarly, a tight proportionality between the critical and the flashing 
temperatures is observed in the right panel of Fig. HI A least-squares fit to 
a single slope parameter results in the relation Tc = 1.28T/, which describes 
well the low temperature points. Note that this constant of proportionality 
is close to the Jaqaman model prediction, Eq. ( 12T1) . and suggests that the 
relation ^ = 2^ is fulfilled. The somewhat poor description of the high 
temperature points can be substantially improved by allowing for an offset. 
The linear relation = —5.39 + l.TlTj gives an excellent description of all 
the data. In this case, the slope is closer to that predicted by the original 
Kapusta model. Note, however, that the 5 highest critical temperatures 
correspond to forces with m* = m or with a poor description of asymmetric 



matter according to Ref. [30|. Modern Skyrme and Gogny forces seem to 



favor the regions ~ 14 - 17 MeV and ~ 11 - 13 MeV. 
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The tight hnear correlation between and Tf provides a clean connec- 
tion between two seemingly unconnected points of the liquid-gas transition 
phase diagram. The correlation is generic, in the sense that it is valid for all 
the mean-fields considered. These findings exemplify the potential of using 
comprehensive sets of interactions to assess the physically allowed param- 
eter spaces of observables. The correspondence between Tf and might 
be useful in relating flashing and critical temperatures in other many-body 
approaches, particularly if one of the two temperatures cannot be accessed 



37( 1 . Note that theoretical correlations between limiting temperatures 



of nuclei and critical temperatures in nuclear matter have also been found 



18l . |31| . These differ from the correlations presented here in that they are 
found after asymmetry, surface and Coulomb effects are taken into account. 

A more detailed analysis of the self-consistent mean-field data is presented 
in Table (H where the saturation and critical properties of the effective forces 
considered are presented in decreasing order of their critical temperature. 
The largest critical temperatures correspond to forces which have at the 
same time a relatively large compressibility and ^ ~ 1. Both inputs are 
therefore potential candidates for physically relevant quantities in a reliable 
determination of the critical temperature. Apart from this, it is hard to find 
any correlation in the data presented in the table. The critical compressibility 
factor, 7c, is relatively stable to the change of mean-field, generally of the or- 
der 7c ~ 0.28 — 0.3. These values are well below the van der Waals prediction, 
but in agreement with experimental data obtained from multi-fragmentation 



58[ and with more sophisticated many-body approaches |37|, |38 



reactions 

Gogny forces are shown separately in the lower end of the Table. In- 
teractions with apparently very similar saturation properties, such as D250 
and DIP, lead to critical temperatures with large differences. This can be 
taken as an indication that the difference in the momentum dependence [i.e. 
in their effective masses) is relevant for critical liquid-gas phenomena. It is 
also interesting to note, by comparing sets D250 to D300, that there is no 
monotonous dependence with the compressibility. 

A comparison between these critical temperatures and those predicted 
by the models presented in the previous Section should be helpful in estab- 
lishing which are the relevant ingredients in the determination of the critical 
properties. A numerical experiment has been performed in which the pre- 
dictions of the critical temperature for each model have been obtained from 
the saturation properties of the Skyrme and Gogny forces and compared to 
the Tc obtained with the self-consistent mean-field data. The results of these 
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experiments are summarized in Figs. M and [TJ which show the "real" 
against the predictions of the different models. Good predictions should line 
up close to the diagonal, which is highlighted with a dotted line. 
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Table 1: Saturation (columns 2 to 5) and critical prop- 
erties (columns 6 to 8) for different Skyrme and Gogny 
interactions in decreasing order of T^. 



Fig. E] focuses on the difference between the critical temperature obtained 
with the Jaqaman model, Eq. f|T8l) . and that of the self-consistent calcula- 
tions. A direct comparison with the pressure of Eq. (1141) can only be estab- 
lished for Skyrme mean-fields. More specifically, in Fig. [5] the critical point 
obtained in Skyrme self-consistent calculations has been compared to those 
obtained by using the Jaqaman predictions with the corresponding values of 
tojt^ and a. The Jaqaman critical temperature is always larger than that of 
the self-consistent mean-field model. For effective interactions with m* = m, 
the self-consistent follows the trend of the Jaqaman model critical tem- 
perature, but it is always shifted upwards by about 5 — 6 MeV. Note that, for 
these interactions, the pressure is not influenced by the momentum depen- 
dent terms of the mean-field and, in the semiclassical limit, will presumably 
be well described by Eq. f|T^ . For the mean-fields with momentum depen- 
dence, however, the discrepancies are quite larger. The closest prediction 
is 8 MeV off, and some differences are as large as 80 MeV. Such a large 
discrepancy indicates that this model is unable to reproduce the thermody- 
namical properties of the self-consistent case. Extensions of this model to 
include effective mass effects as well as degeneracy corrections were already 
discussed in Ref. 5J| and might presumably improve the agreement between 
the different predictions of Tc. 

The results for the original Kapusta model are shown in the left panel 
of Fig. El The predictions for the subgroup of Skyrme forces with m* = m 
(squares) are extremely good, falling very close to the diagonal. For the 
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Figure 5: Self-consistent Skyrme mean- field critical temperature versus the critical tem- 
perature predicted by the Jaqaman model. 



remaining subsets of Skyrme and Gogny forces, Eq. fl26|) produces substan- 
tially worse results and overestimates the critical temperature by up to 9 
MeV. The Kapusta prediction for the critical temperature is therefore valid 
for mean-fields which have an almost quadratic and temperature independent 
momentum dependence. For Skyrme forces, the terms proportional to ^1,^2 
and X2 determine, among other things, the effective mass and are presumably 
responsible for the worsening observed when mg 7^ m. 

As discussed in the previous Section, the effective mass is accounted for in 
a very basic way in the Kapusta model. The differences between the different 
subgroups of Skyrme forces, with and without effective masses, highlight its 
importance in hot nuclear systems. In particular, the effective mass has 
a double effect on the thermodynamical properties. On the one hand, it 
influences their density dependence, as can be seen by comparing Eqs. fl22l) 
and fl52]) . On the other, for a fixed density, the effective mass also affects the 
temperature dependence, i.e. the variations of thermodynamical properties 
with temperature are different for different effective masses. Since the critical 
properties are determined by the density derivatives and the temperature 
dependence of the equation of state, they are particularly sensitive to the 
density dependence of m*. 

The Lattimer-Swesty and Natowitz predictions, Eqs. (|5U]) and ([SI]), have 
been obtained from fits to sets of theoretical predictions to the critical tem- 
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Figure 6: Self-consistent mean-field critical temperature versus the critical temperature 
predictions of the Kapusta (left panel), Lattimer-Swesty (central panel) and Natowitz 
(right panel) models. 



perature. Their accuracy will therefore be limited by the original data used 
to determine the coefficients Cls and Cn- One might, however, try to push 
these parameterizations beyond their original purpose and compare their 
predictions to the mean-field critical points. The results of such analysis are 
presented in the central and right panels of Fig. [61 On average, it appears 
that the Lattimer-Swesty predictions are better (closer to the diagonal), than 
the Natowitz predictions are. For any given interaction, the difference be- 
tween the two predictions is caused by (a) the difference in the coefficients 
and (b) the effective mass at saturation, rriQ. For the Skyrme forces with 
m* = m, the latter factor is not present and, since Cn < Cls, the pre- 
dictions for the critical temperature of the Natowitz model underestimate 
the critical temperature more than the Lattimer-Swesty predictions do. For 
the Skyrme mean-fields with effective masses, the Lattimer-Swesty and the 
Natowitz predictions represent an improvement with respect to the Kapusta 
model. All the Lattimer-Swesty results fall below the diagonal, but they 
overestimate the critical temperature only by about 2 MeV on average (to 
be compared to almost 5 MeV for the Kapusta model). The Natowitz for- 
mula yields predictions which are equally distributed at either side of the 
diagonal, but with a slightly larger deviation from the central results. By 
looking at these figures, it is hard to decide which of the two models produces 
a better agreement with the mean-field predictions. 

The extended Kapusta model yields a qualitatively better description of 
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Figure 7: Self-consistent mean-field critical temperature versus the critical temperature 
predictions of the extended Kapusta model. The inset focuses on the region of low critical 
temperatures. 



the critical temperature in terms of saturation properties. In Fig. U\ the 
results of Eq. (|39i) are compared to mean-field critical temperatures. For all 
the Skyrme forces considered, the agreement between the predicted and the 
mean- field is almost perfect, good to less than 1 MeV accuracy. Note 
that for the Skyrme forces with m* = m, the extended model predictions are 
the same as those predicted by the original Kapusta model. For the forces 
that have a density- dependent effective mass, the extended model supposes 
a substantial improvement with respect to any of the previously discussed 
models. For the forces with > 16 MeV, in particular, the predictions 
follow closely the Skyrme mean-field results. Below this critical temperature, 
the extended Kapusta results show a slight tendency to underestimate the 
critical temperature (see inset), but the results remain always close to the 
diagonal. Comparing these results to other models, it is clear that the density 
dependence of the effective mass is a capital ingredient in the determination 
of the critical temperature. 

The predictions of the extended Kapusta model are worse for the Gogny 
interaction. The maximum observed deviation is of ~ 4 MeV, much larger 
than the largest deviation observed for the Skyrme mean-fields. These dif- 
ferences, however, are similar to those of the Lattimer-Swesty and Natowitz 
predictions for Gogny forces. The sources of the discrepancy might be at- 
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tributed to several factors. To begin with, the density dependence of the 
effective mass for Gogny interactions is not given by that of Eq. More- 
over, the more reahstic momentum dependence introduce modifications in the 
temperature dependence of the pressure that might not be captured properly 
in the model. These considerations challenge the validity of Eq. f l32|) . On 
this account, the extended Kapusta model might need to be reconsidered for 
mean-fields other than those of the Skyrme type. 



In line with the findings of Ref . [3l| , the agreement between the extended 



Kapusta model and the mean-field data suggests that an accurate determi- 
nation of nuclear limiting temperatures might be used to derive saturation 
properties of nuclear matter and, perhaps, identify the density dependence 
of the effective mass. No matter what theoretical framework is chosen to de- 
scribe the phase transition, nuclear mean-field dependences should be pinned 
down carefully. Universal behaviors, such as the connection between Tf and 
Tc observed here, should be identified with a comprehensive set of effective 
interactions. A proper treatment of shell, deformation and continuum con- 
tributions at finite temperature also seems to be a necessary ingredient for 
successful descriptions of data. More importantly, a treatment of clusteriza- 
tion and thermal fluctuations beyond the oversimplified mean-field picture is 
necessary to provide a realistic description of nuclear matter close to a phase 
transition. 



5. Critical exponents 

Critical exponents characterize the properties of phase transitions in a 
universal way [H, [i], |5^, 60 1. In the nuclear physics context, they have been 



introduced to analyze yields of multi-fragmentation reactions and the onset 
of percolation jHs], [HH. The critical exponents of infinite nuclear matter have 
received far less attention. Calculations can easily be performed in simplified 
models, as those presented in Section [31 but they are harder to implement 
numerically for self-consistent mean-field models. Some critical exponents 
have been computed for a single Skyrme mean-field in Ref. [62| to address 
issues related to bimodality in isospin asymmetric matter. 

The calculations of critical exponents for different mean-fields are of in- 
terest for theoretical reasons. On the one hand, simplified models, as has 
been discussed in the two previous sections, might not be able to reproduce 
the liquid-gas phase transition properties of mean-field (or more complicated 
many-body) calculations. Consequently, their predictions for critical expo- 
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nents might not be rehable. Extensive mean-field calculations are needed to 
confirm numerically the values of the exponents and to check its indepen- 
dence from the underlying effective interaction. It is particularly interesting 
to note that the self-consistency at the mean-field level can only be imposed 
numerically and thus precludes any analytical development. On the other 
hand, it is well-known from magnetic systems that critical exponents are 



sensitive to the range of the interaction [59|. Therefore, it is plausible that 
self-consistent mean-field calculations of homogeneous nuclear matter have 
different critical exponents for interactions which have either zero or finite 
range. 

Critical exponents give unique information on the behavior of the ther- 
modynamical properties of nuclear matter close to the critical point. A nu- 
merical code has been implemented to compute three of these exponents 
for mean-field calculations of nuclear matter. The order parameter in the 
liquid-gas phase transition is given, at each temperature, by the difference 
of the liquid and the gas phase densities, i. e. the coexistence line. The first 
exponent, /3, characterizes how the order parameter approaches the critical 
point ^: 

Pi-Pa-^ i-rr , (40) 

where the reduced temperature, r = is a measure of the distance 

to the critical point. The critical exponent 7 is given by the isothermal 
compressibility. 



1 dp 
Kt ^ dp 



(41) 



close to the critical temperatur^, 

Kt^{-t)-^. (42) 



^Strictly speaking, one should consider different critical exponents for the liquid and 
gas branches at r < 0, so that (/3,7, (5) would become {Pi, ^1,61) and il3g,'yg,Sg). The 
numerical results presented below fulfill I3g = f^i, jg — 7; and Sg = Si to within 5%, 
6% and 10% accuracy, respectively. The results that are quoted have been obtained by 
averaging liquid and gas branches to eliminate uncertainties in the determination of pc- 
Moreover, for r > and p = pc one can also define a 7 critical exponent ^pd\. Scaling laws 
suggest 7 = 7', which is fulfilled numerically to a 1% accuracy. 
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Moreover, the density dependence of the critical isotherm is described by the 
critical exponent, 6^: 

\P-Pc\ ~ \p~Pc\'. (43) 

Finally, for nuclear matter, the heat capacity at p = pc also approaches the 
critical point with a power-law behavior: 

^v = T^- ■ (44) 

a will not be explicitly computed here, but its value could be derived from 
the other critical exponents by means of either of the following scaling law 



relations 60 



2-a = /3(5 + l), 2-a = 7 + 2/3. (45) 

The critical exponent /3 is closely related to the temperature dependence 
of the coexistence curve. A basic property of the liquid-gas phase transition 
for atomic substances is the law of corresponding states, which dictates that 
the equation of state of simple atomic liquids and gases is the same after 



normalizing it to the respective critical properties |63[. An important 
consequence of this law is that, after rescaling by the critical density and 
temperature, the coexistence curve of a wide range of substances is the same. 
It is important to note that this is an empirical law, which indicates that the 
thermodynamics of all these substances is similar. In the nuclear matter case, 
one might wonder whether different coexistence lines, obtained from different 
mean-fields, follow a similar law of corresponding states. Analogous studies 
have been carried out in the context of molecular dynamics simulations of 



classical liquids |42|. For that particular case, the coexistence curve obtained 
with different inter-atomic potentials fulfills a law of corresponding states. 

The scaled coexistence curves for a few representative mean-fields are 
shown in Fig. [HI SV (long-dash-dotted line) and SLy9 (short-dash-dotted 
line) represent extreme behaviors of the coexistence curves. The remaining 
mean-field coexistence curves fall within these two results. While in the gas 
phase there seems to be a certain degree of agreement between different mean- 
fields, substantial differences are found in the liquid branch even for values 
of T/Tc ~ 0.9, relatively close to the critical point. The breakdown of the 
law of corresponding states provides an interesting insight on the interaction 
dependence of the thermodynamical properties. To begin with, it indicates 
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Coexistence curve 




Figure 8: Normalized coexistence curve for different mean-fields. The small circles corre- 
spond to the Guggenheim parameterization, Eq. (j46|) . The diameters of each mean- field 
are also shown. 



that the equation of state of different mean-fields is not simply related to one 
another by a simple scaling with the critical properties. This is in contrast 
to classical liquids, in which plausibly the differences in interactions are less 



acute and such a scaling is accurate 42|. As a matter of fact, on physical 



grounds, it is only the thermodynamics close to the critical point that should 
be interaction independent and dictated by the critical exponents. The mean- 
field predictions of Fig. |8] are compared to the Guggenheim parameterization: 

Pc 4 4^^' ^ ^ 



where the -|- (— ) sign is used for the liquid (gas) branch 63|]. Instead of the 



original /3 = |, which is suitable for classical gases |60|, the value (3 = ^ has 



been used, in accordance to the results discussed in the following paragraphs. 
The Guggenheim formula produces a somewhat average behavior, falling 
within all the mean-field results and might be used for phenomenological 
applications. 

Useful information regarding the phase coexistence can also be gathered 
by analyzing the so-called "diameter" of the transition, i.e. the average of 
the coexistence densities at each temperature. The diameter is given by the 
almost vertical, central line shown in Fig. El The tilt of the diameter with 
respect to the vertical indicates that the coexistence curve is not perfectly 
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Figure 9: Isothermal compressibility as a function of reduced temperature for different 
mean-fields. Below Tc, the liquid and gas branches are obtained at the corresponding 
coexistence densities. Above Tc, Kt is computed at p^. 



symmetric. As a matter of fact, the law of rectilinear diameter sates that the 
diameter scales linearly with r: 



where c is some numerical coefficient [42|, |60|, 163|. As observed in Fig. |8l 
the diameter is indeed rectilinear for all the mean-fields considered within a 
wide range of temperatures (the linear behavior breaks down very close to 
Tc, with a deviation related to a [i^]). Every mean- field predicts a different 
tilt {i.e. a different value of c) of the diameter. Some of these agree with 
the value c = |, predicted by Eq. fH^ . but generally lower values (as low as 
c ~ 0.13) are found. These indicate that mean- field coexistence curves are 
somewhat more symmetric than the Guggenheim parameterization. 

The isothermal compressibility diverges close to the critical point with 
a power law behavior governed by the 7 critical exponent. Kq^ is shown in 
Fig. [9] as a function of the reduced temperature for different mean-fields. 
Below Tc, the gas and the liquid branch are shown independently. A fairly 
large variation of is observed for the different mean-fields in the gas 
region (note the logarithmic scale), in contrast to the relatively similar gas 
coexistence curves seen in Fig. [HI For the liquid branch, the predictions for 
Kt^ are relatively closer to each other. It is important to note that i^y is not 
related to the nuclear compressibility, A'q- Only in the T — limit, the value 
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Figure 10: Logarithmic plots used to determine the critical exponents. Left upper panel: 
difference between liquid and gas coexistence densities versus reduced temperatures. Left 
lower panel: isothermal compressibility versus reduced temperatures. Right panel: critical 
isotherm versus density. Symbols correspond to mean-field data, while lines are the results 
of linear regressions. 



of Kt in the liquid branch is given by the compressibility and the saturation 
density, Kt = Forces with larger compressibilities, such as ZRl and 

SV, will therefore lead generally to lower Kt^s. Above the critical point, Kt 
is taken at p = pc and its value decreases steadily, in accordance to the Kt ~ 
behavior. Again, a certain degree of spread between different mean- 
field predictions is found, particularly among those with different nuclear 
compressibilities. Note that, in spite of these differences, the critical exponent 
7 is the same for all the interactions considered (see below). 

For systems with no thermal fluctuations, the Landau mean-field theory 
of phase transition predicts the following values for the critical exponents 
3: 

/3 = ^, 7 = 1, S = 3, a = 0. (48) 
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A first indication of the nuclear matter critical exponents can be obtained by 
the models introduced in Section El Indeed, the Jaqaman model reproduces 



the Landau theory results [5J, |65[ and one can easily show that the Kapusta 
model (and its extension) also fall in the same category. As mentioned earlier, 
however, these basic parameterizations of the equation of state are too simple 
to describe the liquid-gas coexistence region of nuclear matter derived from 
self-consistent calculations. The latter provide a consistent description of hot 
nuclear matter, with no assumptions whatsoever on temperature or density 
dependences, providing a more sound analysis from a theoretical perspective. 
It is clear, however, that the mean-field description has important limitations 
in describing nuclear liquid-gas phenomena. As a consequence, the following 
results concerning critical exponents should not be taken as realistic, but 
rather as a theoretical confirmation of the validity of Landau mean-field 
theory for homogeneous nuclear matter. 

The self-consistent calculations give critical exponents of the Landau 
mean-field type. As expected from physical grounds, this result is inde- 
pendent of the effective interaction under consideration. To prove the latter 
point. Fig. [TUl shows the logarithmic plots necessary to determine the values of 
/3, 7 and 6. The data displays a perfectly linear behavior. The slopes of these 
lines correspond to the critical exponents, and they are indeed independent 
of the interaction. In agreement with the previous considerations, however, 
the offsets of the lines do depend to a certain extent on the interaction. 

It is important to indicate that, for an accurate calculation of these expo- 
nents, the coexistence line needs to be determined reliably for temperatures 
very close to the critical point. In addition, both the critical temperature 
and density need to be known with relative precisions of at least 4 significant 
digits. While this is of a lesser concern for Skyrme forces, the calculations for 
finite-range Gogny mean-fields have proven to be harder and involve longer 
computational times. For these interactions, coexistence curves have been 
determined typically down to r ~ 10~^ and, as a consequence, somewhat 
less accurate results have been obtained. Note that the data for Kt has been 
obtained as an average of the gas and the liquid branches. Similarly, for the 
critical isotherm, densities above and below pc have been averaged. This pro- 
cedure cancels out part of the numerical uncertainties in the determination 
of the critical density and temperature. 

The critical exponents have been obtained numerically by a linear regres- 
sion of the mean-field data presented in Fig. [TOl To determine /3 and 7, a 
logarithmic mesh of 12 values between r = 10~^'^ and 10~^'^ (10~^) has been 
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/3 


7 


S 


ZRl 


0.4996±0.0001 


0.9951±0.0013 


3.0000±0.0001 


DIM 


0.5019±0.0017 


1.0074±0.0103 


3.0124±0.0034 


BSkir 


0.4993±0.0002 


0.9934±0.0017 


3.0006±0.0002 


SLy9 


0.4989±0.0003 


0.9926±0.0020 


3.0007±0.0002 


SV 


0.5001±0.0006 


0.9978±0.0036 


3.0067±0.0005 



Table 2: Critical exponents obtained from linear regression fits of the results of Fig. [TOl 
Quoted errors correspond to the standard deviations of the slope. 



used for Skyrme (Gogny) forces. Similarly, to find 6, 12 logarithmic points 
have been distributed between = 10~°'^ and 10~^'^^. The results of the 
slopes of the corresponding regressions for the five interactions considered in 
this Section are presented in Table |2l The quoted errors correspond to the 
standard deviations of the slope and can be taken as a measure of the quality 
of the fits. The extremely small values of these errors provide an indication 
of the very good linear behavior of the data in the logarithmic scales and 
therefore validate the power law behavior of the self-consistent mean-field 
data. The agreement with the Landau predictions is excellent, accurate to 
the third digit in most cases. Minute deviations from this behavior should 
plausibly be attributed to numerical errors associated with the limited accu- 
racy of the self-consistent numerics close to the critical point. While only 5 
mean-fields are shown here for practical purposes, identical results have been 
obtained for the remaining 36 Skyrme and 7 Gogny interactions considered 
in the previous Section. 

For real physical systems, the values of the critical exponents can devi- 
ate substantially from the Landau mean-field predictions j2|, |6o[. Yields in 
multifragmentation reactions, for instance, can be described with /3 ~ 0.3 



and 7 ~ 1.4 [66|. The fact that Landau-type exponents are found for 
Hartree-Fock calculations of nuclear matter is a sign of its inability to prop- 
erly describe the critical properties. The self-consistent mean-field approach 
represents a minimization of the thermodynamical potential, subject to the 
constraint that the density matrix is of the one-body type 23|]. The physi- 
cal trajectories that this minimization procedure yields are unique and have 
small thermal fluctuations by construction. Critical exponents different from 
Landau theory, however, can only be achieved if thermal fluctuations, pre- 
sumably larger than those introduced by one-body approximations, are taken 
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into account 



It is plausible that this can be achieved by either consid- 
ering correlated many-body approaches or by including additional stochastic 
fluctuations on top of the mean-field 671 . 

To summarize, the normalized equations of state, coexistence curves and 
thermodynamical properties of nuclear matter coming from different self- 
consistent mean-field calculations do not follow a law of corresponding states. 
At low temperatures, the thermodynamical properties are actually sensitive 
to details of the interaction, their ranges and structures. In spite of these 
differences, critical exponents are the same, independently of the mean-field. 
The physics close to the critical point is therefore not affected by the structure 
of the short-range nuclear effective interaction and, at the critical point itself, 
the long-range thermal fluctuations dominate the thermodynamics of nuclear 
matter. For self-consistent mean-field calculations of infinite matter, the 
numerical values of the critical exponents agree with those given by Landau 
mean-field theory. 



6. Conclusions and outline 

Hot nuclear matter has been studied within the self-consistent Hartree- 
Fock approach using different mean-field interactions of the Skyrme and the 
Gogny types. The basic aim of this work has been to provide a systematic 
analysis of the liquid-gas phase transition in the mean-field approximation for 
nuclear matter. Calculations performed with several mean- fields are helpful 
in elucidating which are the physically relevant properties for the liquid-gas 
transition, since they provide means of eliminating the mean-field dependence 
of the results. Subtracting this dependence, a more systematic understand- 
ing of the properties of the liquid-gas phase transition is obtained. This 
procedure might also be necessary in more realistic descriptions of nuclear 
critical phenomena. To my knowledge, this is the first time that a generic 
analysis with different effective interactions within a self-consistent theory is 
used to this end. 

A step-by-step improvement of our understanding of the temperature 
dependence of nuclear properties requires a full control of symmetric nu- 
clear matter as a first, initial step. Two points stand out in the phase di- 
agram due to their physical relevance: the saturation point and the critical 
point. Liquid-gas critical points have been computed for different Skyrme 
and Gogny interactions and correlations with zero temperature saturation 
properties have been investigated. Parameter spaces for physically allowed 
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critical densities and temperatures have been identified by using a compre- 
hensive set of mean-field calculations. On the one hand, it is difficult to find 
clear connections between saturation and finite temperature properties. On 
the other hand, a tight linear relation has been found between the flashing 
and the critical density, as well as between the flashing and critical tem- 
peratures. Given the physical equivalence between the flashing temperature 
of nuclear matter and the limiting temperature of nuclei, this correlation 
might be useful in relating the nuclear matter critical point to the flashing 
properties of nuclei found in reactions. 

A qualitative understanding of the properties of the phase transition can 
be obtained by using analytical parameterizations of the equation of state. 
The Jaqaman and the Kapusta models represent two of these parameteri- 
zations and arc among the most widely used in the literature. The critical 
point can be solved explicitly in these models, providing relations between 
the critical point and parameters related to the zero-temperature equation 
of state. These models also yield correlations between the flashing and the 
critical points, in overall agreement with those found with the mean-field re- 
sults. The temperature dependence of the Jaqaman and the Kapusta models 
is based on degeneracy expansions, either in the semi-classical or the de- 
generate regimes. Unfortunately, the liquid-gas phase transition in nuclear 
matter occurs between these two regimes and none of these models is able 
to predict reliably the mean- field results of the critical point. In particular, 
the critical temperatures are generally overestimated by few to some tenths 
of MeV. Lattimer-Swesty and Natowitz et al. have proposed phenomenolog- 
ical modifications to improve the predictive power of the Kapusta model. 
Both formulae do a better job in predicting the critical temperature of the 
Skyrme and Gogny forces, but fail to give quantitatively accurate results 
when compared to mean-field systematics. 

A common feature of the analytical models is a very crude treatment of 
the effective mass, in spite of the fact that m* has been identified as a basic 
ingredient in the physics of hot nuclei and the liquid-gas phase transition. 
The Kapusta model has been extended to include the density dependence of 
the Skyrme effective mass. While this model cannot be solved analytically, an 
expansion close to the critical point yields systematic corrections to the Ka- 
pusta model predictions. At first order, while the critical density is basically 
unaffected, the critical temperature is modified by a multiplicative factor 
that is equivalent to the effective mass at a suitable sub-saturation density. 
This expression provides an excellent prediction of the critical temperature 
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for all the Skyrme mean-fields, to within 1 MeV accuracy, and represents a 
quantitative improvement with respect to the other models. These results 
illustrate the importance of the density dependence of the effective mass, 
i.e. the momentum dependence, in mean-field models at finite temperature. 
The extended model works worse for Gogny forces, with a tendency to un- 
derestimate the critical temperature by up to 4 MeV. It is plausible that a 
modification is needed to account for the different density dependence of the 
effective mass with this mean-field. 

As an additional outcome to this study, the (3, 7 and 6 critical exponents 
of the transition have been computed numerically. As expected from general 
physical grounds, no interaction dependence is found in the critical exponents 
which, even for the finite-range Gogny interaction, agree with the predictions 
of the Landau mean-field theory of critical phenomena to less than a 1% ac- 
curacy. The lack of effective interaction dependence of the critical exponents 
is in contrast to the thermodynamical properties, which do depend on the 
mean- field interaction. As a matter of fact, if the coexistence curve of dif- 
ferent mean-fields is rescaled to the corresponding critical properties, a clear 
disagreement is found i.e. the law of corresponding states is not valid for 
different nuclear mean-fields. 

Further generalizations of this work to asymmetric nuclear matter and 
finite nuclei would provide a consistent route to connect the effective in- 
medium nuclear interaction and the finite temperature properties of nuclear 
matter. For isospin imbalanced systems, for instance, the critical temper- 
ature might be correlated to the symmetry energy (and/or its derivatives). 
Studies with comprehensive mean-field sets, analogous to the analysis pre- 
sented here, could easily test this idea. Similarly, one can nowadays carry 
out fully self-consistent mean-field calculations of semi-infinite matter at fi- 



nite temperature [68|. From these studies, one could determine the tem- 
perature and isospin dependence of nuclear surface properties. At relatively 
low temperatures, moreover, it could be interesting to analyze the temper- 
ature dependence of nuclear structural properties. Mean-field (or, for this 
purpose, energy density functional) nuclear structure calculations have pro- 
vided and continue to provide accurate predictions for single-particle and 
bulk properties in wide regions of the nuclear chart. Their extension to finite 
temperature yields a theoretically consistent description of hot nuclei that 
could, for instance, yield predictions for the isotopic dependence of limiting 
temperatures. A large source of uncertainty in such analysis, however, is 
the extremely important role played by Coulomb forces. For finite systems. 
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this is well beyond the reach of mean-field theory and additional theoreti- 
cal developments might be needed to provide a realistic description of the 
coexistence and spinodal regions. 
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